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We study a gravitational model with curvature-squared R? and curvature-quartic R 4 nonlineari- 
ties. The effective scalar degree of freedom <j> (scalaron) has a multi- valued potential U((j>) consisting 
of a number of branches. These branches are fitted with each other in the branching and monotonic 
points. In the case of four-dimensional space-time, we show that the monotonic points are pene- 
trable for scalaron while in the vicinity of the branching points scalaron has the bouncing behavior 
and cannot cross these points. Moreover, there are branching points where scalaron bounces an 
infinite number of times with decreasing amplitude and the Universe asymptotically approaches 
the de Sitter stage. Such accelerating behavior we call bouncing inflation. For this accelerating 
expansion there is no need for original potential U ((f>) to have a minimum or to check the slow-roll 
conditions. A necessary condition for such inflation is the existence of the branching points. This 
is a new type of inflation. We show that bouncing inflation takes place both in the Einstein and 
Brans-Dicke frames. 
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I. INTRODUCTION 

Starting from the pioneering paper the nonlin- 
ear (with respect to the scalar curvature R) theories 
of gravity f(R) have attracted the great deal of inter- 
est because these models can provide a natural mech- 
anism of the early inflation. Nonlinear models may 
arise either due to quantum fluctuations of matter 
fields including gravity [2| , or as a result of compactifi- 
cation of extra spatial dimensions Q. Compared, e.g., 
to others higher-order gravity theories, f(R) theories 
are free of ghosts and of Ostrogradski instabilities [J] . 
Recently, it was realized that these models can also ex- 
plain the late time acceleration of the Universe. This 
fact resulted in a new wave of papers devoted to this 
topic (see e.g., recent reviews 

The most simple, and, consequently, the most stud- 
ied models are polynomials of R: f(R) — Yl n =o C n R n 
(k > 1), e.g., quadratic R + R 2 and quartic R + R 4 
ones. Active investigation of these models, which 
started in 80-thyears of the last century 0, [8| , contin- 
ues up to now . Obviously, the correction terms (to 
the Einstein action) with n > 1 give the main contri- 
bution in the case of large R, e.g., in the early stages 
of the Universe evolution. As it was shown first in [1| 
for the quadratic model, such modification of gravity 
results in early inflation. From the other hand, func- 
tion f(R) may also contain negative degrees of R. For 
example, the simplest model is R + R^ 1 . In this case 
the correction term plays the main role for small R, 
e.g., at the late stage of the Universe evolution (see 
e.g. [13, HU and numerous references therein). Such 
modification of gravity may result in the late-time ac- 
celeration of our Universe [12| • Nonlinear models with 
polynomial as well as i? _1 -type correction terms have 
also been generalized to the multidimensional case 
(see e.g., Qjl ED, EMI). 

It is well known that nonlinear models are equiva- 
lent to linear-curvature models with additional scalar 



field <p (dubbed scalaron in [Tj]). This scalar field cor- 
responds to additional degree of freedom of nonlinear 
models. The dynamics of this field (as well as a pos- 
sibility of inflation of the Universe) is defined by po- 
tential U{R{4>),4>) (see Eq. gp) below 1 ) where R = 
R{4>) is a solution of Eq. flOf : exp (A<f>) = df/dR. 
Usually, models or a particular cases of these models 
are considered where this equation has only one solu- 
tion. In this case, potential U is a one- valued function 
of <J). However, in the most general case this equa- 
tion has more than one solution and potential becomes 
a multi-valued function with a number of branching 
points (see e.g., [HI])- Investigation of the dynamical 
behavior of scalar field and Universe in such models 
(especially in the vicinity of the branching points) is 
not a trivial problem and may result in new impor- 
tant effects. Therefore, it is of interest to consider the 
models with multi- valued potentials. 

In the present paper, we study an example of such 
models. Here, f(R) has both quadratic R 2 and quar- 
tic i? 4 contributions. In this case Eq. (|2.3|) is a cubic 
equation with respect to R and may have, in general, 
three real solutions/branches Ri((f>) (i = 1. 2, 3). We 
have investigated this model in our paper [16| . How- 
ever, in this paper we considered a special case of one 
real solution in D = 8 space-time. Now, we study the 
most interesting case of three real solutions. These 
solutions are fitted with each other in the branch- 
ing points. There are also another type of matching 
points where one-valued solutions are fitted with the 
three- valued solutions. In the vicinity of these points 
potential U(<j>) is a monotonic function. Thus, these 
latter points dubbed monotonic ones. The main aim 
of the paper consists in the investigation of dynamical 
behavior of the system in four dimensional space-time 
in the vicinity of the branching and monotonic points. 
We show that dynamics is quite different for branch- 
ing and monotonic points. The monotonic points are 
penetrable for scalaron while in the branching points 
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Starting from Sec. II, we denote the scalar curvature of the 
original nonlinear model by R. 
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scalaron has bouncing behavior. There are branch- 
ing points where scalaron bounces an infinite number 
of times with decreasing amplitude and the Universe 
asymptotically approaches the de Sitter stage. Such 
accelerating behavior we call bouncing inflation. We 
should note that for this type of inflation there is no 
need for original potential U(<p) to have a minimum 
or to check the slow-roll conditions. A necessary con- 
dition for such inflation is the existence of the branch- 
ing points. This is a new type of inflation. We show 
that this inflation takes place both in the Einstein 
and Brans-Dicke frames. This is the main result of 
our paper. We think that scalaron field and the Uni- 
verse have the similar behavior in the vicinity of the 
branching points for others polynomial and i? _1 -type 
models resulting in both early inflation and late-time 
acceleration. Of course, it is necessary to conduct ad- 
ditional studies of these models, to confirm or refute 
this assertion. 

The paper is structured as follows. In Sec. II 
we study briefly the equivalence between an arbi- 
trary nonlinear f(R) theory and theory linear in an- 
other scalar curvature R but which contains a scalaron 
field cf>. In Sec. Ill we consider a particular exam- 
ple of nonlinear model with curvature-quadratic and 
curvature-quartic correction terms and obtain solu- 
tions/branches Ri((f>). The fitting procedure for these 
branches is proposed in Sec. IV. The dynamics of the 
scalaron and the Universe is investigated in Sec. V. 
Here, we parameterize the scalaron potential in such 
a way that it becomes a one-valued function. It gives 
a possibility to study the dynamical behavior of the 
system in the vicinity of the branching and monotonic 
points. We show that in the vicinity of the branching 
point the scalaron field bounces an infinite number of 
times with decreasing amplitude and the Universe ac- 
quires the accelerating expansion approaching asymp- 
totically to the de Sitter stage. Such accelerating ex- 
pansion we call bouncing inflation. A brief discussion 
of the obtained results is presented in the concluding 
Sec. VI. In Appendix A, we show that bouncing infla- 
tion in the vicinity of the branching point takes place 
also in the Brans-Dicke frame. 



where 



dR 



> 



.4 



D 



D - 1 



(2.3) 



and where the self- interaction potential £/(</>) of the 
scalar field <fr is given by 



U{4>) 
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(fT D/{D ~ 2) [Rf'-f] , (2-4) 
e- B *[R^)e A *-f{R{<f))] , (2.5) 

(2.6) 



D 



y/(D-2)(D-l) 



The metrics g a b, gab and the sc alar curvatures R, R 
of the two theories (|2.1I) and (|2.2[) are conformally 
connected by the relations 2 
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R 



(f) 2/{2 - D) [R+^(fr 2 g ab daf'dbf' 
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via the scalar field <f> = ln[/'(i?)]/A This scalar field 
(f>, known as scalaron carries an additional degree 
of freedom of original nonlinear model. 

According to our definition (|2.3I) . we consider the 
positive branch f'(R) > 0. Although the negative 
/' < branch can be considered as well (see e.g. Refs. 
H, E3, EH)- However, negative values of f'(R) result 
in negative effective gravitational " constant" G e / / = 
K 2 D /f. Thus /' should be positive for the graviton to 
carry positive kinetic energy (see e.g., [!(). 

From action (|2.2p we obtain the equation of motion 
of scalaron field 6: 



(2.9) 



II. GENERAL SETUP 



If scalaron potential U(<p) has a minimum in a point 
^0 : 



It is well known that nonlinear theories 
1 



S 



% K D JM 



l\f(R) , 



(2.1) 



where f(R) is an arbitrary smooth function of a 
scalar curvature R = R[g] constructed from the 
D— dimensional metric g a b (a, b = 1, . . . , D) are equiv- 
alent to theories which are linear in another scalar 
curvature R but which contains an additional self- 
interactingscalar field. According to standard tech- 
niques 0,H|, the corresponding R— linear theory has 
the action functional: 



dU 



A 



2(Z3 - 2) 
h := Df - 2Rf, 



{n -D,(D-2) h 



= o, 



<Po 



ft(0o)=O (2.10) 



then we can define the mass squared of the scalaron 
0: 



HI,! 



d 2 U 



1 



(D - 2)f - 2Rf" 



^ 2{D-l)f" (f>)V(D-2) 
(D - 2)f - Dff'/f 



2(D- l)/"(/')2/(D-2) 



> 0. 



4>o 
(2.11) 



4>o 



S 



1 



2^£) J M 



d D x^\g~\ [R[g] - g ah ^, a 4>.b 



217(0)] , 
(2-2) 



2 The metrics g a i, and g a b represent the Einstein and Brans- 
Dicke frames, respectively. 
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Similar expression for the scalaron mass squared is 
also given e.g., in [6]. The only difference consists 
in additional conformal prefactor l/(/') 2 /(-°~ 2 ) orig- 
inated from the conformal metric transformation 3 
( 2.7|). Up to this prefactor, the positiveness condition 
(|2.1ip of the mass squared coincides with the stabil- 
ity condition of de Sitter space in f(R) gravity with 
respect to inhomogeneous and homogeneous pertur- 
bations d, I2lj . Additionally, to avoid the Dolgov- 
Kawasaki instability [22| (instability with respect to 
local perturbations), it is also required that f"(R) > 
(see also [23|). 

For further research is useful to introduce a new 
variable 




X := - 1 . 



(2.12) 



Then, potential (|2.6I) and its first derivative read, cor- 
respondingly: 



1 



U(X) = - (X + iy B/A [R(X + 1) - /] (2.13) 



and 



+ \{X + \y B/A R. (2.14) 



To conclude this section, we would like to remind 
that in multidimensional case, to avoid the effective 
four-dimensional fundamental constant variation, it 
is necessary to provide the mechanism of the inter- 
nal spaces stabilization. In these models, the scale 
factors of the internal spaces play the role of addi- 
tional scalar fields (geometrical moduli/gravexcitons 
24]). To achieve their stabilization, an effective po- 
tential should have minima with respect to all scalar 
fields (gravexcitons and scalaron). Our previous anal- 
ysis (see e.g. [25]) shows that for a model of the form 
(|2.2p the stabilization is possible only for the case of 
negative minimum of the potential U(cj>). However, it 
is not difficult to realize that it is impossible to freeze 
out the internal spaces in such AdS universe. Indeed, 
in these models scalar fields decrease their amplitude 
of oscillations around a minimum during the stage of 
expansion of the Universe (due to a friction term in 
dynamical equation of the form of (|5.2p below) until 
the Universe reaches its maximum. Then, the Uni- 
verse turns to the stage of contraction and the am- 
plitudes of scalar fields start to increase again. Thus, 
geometrical moduli are not stabilized in such models. 
Therefore, in our present paper we do not investigate 
the problem of the extra dimension stabilization but 
we focus our attention on the dynamics of the scalaron 
field and the Universe in four-dimensional case. 



3 Conformal transformation for mass squared of scalar fields in 
models with conformally related metrics is discussed in [20ll , 



FIG. 1: The s urfa ces /' = as a functions R — R(a,y) 
for the model (13.111. 



III. THE R 2 + -MODEL 

In this section we analyze a model with curvature- 
quadratic and curvature-quartic correction terms of 
the type 



f(R) = R + aR 2 +jR 4 



2A, 



(3.1) 



We start our investigation for an arbitrary number of 
dimensions D but in the most particular examples we 
shall put D = 4 (unless stated otherwise). First of all, 
we define the relation between the scalar curvature R 
and the scalaron field <fi. According to eq. (|2.3|) we 
have: 



1 + 2aR + 4jR 3 



(3.2) 



The definition (|2.3[) /' = exp(A</>) clearly indicates 
that we choose the positive branch /' > 0. For our 
model (|3.1I) . the surfaces /' = as a functions R = 
R(a, 7) are given in Fig. (flj. As it easily follows from 
Eq. (|3.2I) . points where all three values R, 7 and a 
are positive correspond to the region /' > 0. Thus, 
this picture shows that we have one simply connected 
region /' > and two disconnected regions /' < 0. 
Eq. (|3.2p can be rewritten cquivalently in the form: 



= e M 



X = e 



-00 < q> < +00 
-1 < X < +00 . 



(3.3) 



(3.4) 



Eq. p.3p has three solutions i?i,2.3, where one or three 
of them are real-valued. Let 



q:=— r-=—X 

67 ' 87 



The sign of the discriminant 

Q := r 2 + q 3 



(3.5) 



(3.6) 
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defines the number of real solutions: 



Q > 

= o 

Q < 



3^ = Vi, R 1 = R 2 , 
SRt = Vi . (3.7) 



Physical scalar curvatures correspond to real solu- 
tions Ri(X). It is the most convenient to consider 
Ri = Ri(X) as solution family depending on the two 
additional parameters (a, 7): sign(a) = sign{^) ==> 
Q > 0, sign(a) 7^ sign(j) =>■ Q 0. The case 
sign{a) — sign{~f) was considered in our paper [l6| . 
In the present paper we investigate the most interest- 
ing case signia) ^ sign{^) of multi- valued solutions. 
For Q > the single real solution .Ri is given as 



Ri 



,1/2 



1/3 



r-Q 



1/2 



1/3 



zi + z 2 , (3.8) 



where we can define z\ 2 in the form: 



p 2 =r 2 -Q = -q\ 



cosh((9) 



(3.9) 



Taking into account eq. (13.51) . the function X reads 
X{6) = 87 V /Z g 5 cosh(6>) . (3.10) 

The three real solutions R.1,2, 3 (X) for Q < are 
given as 



Ri 
Ri 

R a 



Sl + s 2 , 



-(-1 + n/3>i + i(-l - iV3)s 2 



'si+e 1 3 s 2 , 



^(-l-n/3>i + ^(-l + h/3) 



•S2 



3 si + e' 3 s 2 , 



(3.11) 



where we can fix the Riemann sheet of Q 1 ^ 2 by setting 
in the definitions of si 2 



Si, 2 



±i\Q\ 



1/2 



1/3 



(3.12) 



A simple Mathematica calculation gives for Vieta's 
relations from p. lip 



R!+R 2 + R 3 = 0, 
R 1 R 2 + R1R3 + R2R3 = -3s lS2 = 3q, 



RiR 2 R 3 



2r. (3.13) 



In order to work with explicitly real-valued Ri we 
rewrite Si i2 from (|3.12l) as follows 



Slj2 = |6|V3 e ±»V3 ; 



H 2 = 



cos(tf) 



\Q\=r - Q = -q A 
r 



\b\ 



-Q 



(3.14) 



and get via p. lip 

: 2|6| 1 / 3 cos(i?/3), (3.15) 
e"^s 2 = 2|6| 1/3 cos (i9/3 + 27t/3), 
^si +e l T S2 = 2|&| 1/3 cos (i?/3-2tt/3) 



-Ri = si + s 2 
#2 = e'^si - 



R 3 = e 
or 



2|6| 



1/3, 



-<7 cos 



2ttA:\ 



3 

■d + 2rrk 



k = -1,0, 1. (3.16) 



In order to understand the qualitative behavior of 
these three real-valued solutions as part of the g loba l 
solution picture we first note that, according to p.3p . 
we may interpret X as single- valued function 



X(R) =4 7 i? 3 + 2aR 



(3.17) 



and look what is happening when we change (a, 7). 
Obviously, the inverse function R(X) has three real- 
valued branches when X(R) is not a monotonic func- 
tion but instead has a minimum and a maximum, i.e. 
when 



d R X := X' = 12 7 i? 2 + 2a = 



R 2 



67 

. (3-18) 

has two real solutions R± = —a/ (67) and corre- 
sponding extrema 4 



X(R ± ) = -aR± . 



(3.19) 



It should hold sign (a) 7^ sign (7) in this case so that 
we find 

7 > 0, a < : X max = X(iL), X mm = X(R+) 
7 < 0, a > : X mQX = X(i? + ), X mm = X{iL.) . 

(3.20) 

The transition from the three-real-solution regime to 
the one-real-solution regime occurs when maximum 
and minimum coalesce at the inflection point 



R+ = R- = 



a = 0, 7 ^ 0. 



(3.21) 



(We note that here we consider the non-degenerate 
case 7^0. Models with 7 = are degenerate ones 
and are characterized by quadratic scalar curvature 
terms only.) Due to — 1 < X < +00 we may consider 
the limit X — > +00 where in leading approximation 



4 7 i? 3 X 



+00 



so that 



R(X — >• 00) — > sign (7) x 00 . 



(3.22) 
(3.23) 



It is worth of noting that f"(R±) = X'(R ± ) = 0. 
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Leaving the restriction X > — 1 for a moment aside, 
we have found that for 017 < there exist three real 
solution branches 7^1,2,3: 

7>0: -oo<Ui<R-, -00 < X < 

R+ < K 3 < +00, X mm <X< +00, 

OO > X > X m i n , 



— X maXl 



>x>x 



rain j 



7 < : -00 <K X <R- 



/?.. < %:i < -+ X.„ 



X — ^-raax 1 



(3.24) 



It remains for each of these branches to check which 
of the solutions Rk from Q3.16P can be fit into this 
scheme. Finally, one will have to set the additional 
restriction X > — 1 on the whole picture. 



IV. THE FITTING PROCEDURE 

We start by considering a concrete example. For 
definiteness, let us assume 7 > 0, a < 0. The pairwise 
fitting of the various solution branches should be per- 
formed at points where Q = and different branches 
of the three-solution sector are fitted with each other 
or to the branches of the one-solution sector. The 
points with Q = correspond to the X— v alues X m i n 
and X max . Explicitly we have from (|3.19[) 



4 4 / a~ 

X(R±) = -aR±=±-aJ-— (4.1) 

and for the concrete configuration 7 > 0, a < 

V • ' 4 



X,, 



X(R + ) = ^aJ^<0. (4.2) 



Next, wc find from the defining equation (|3.14p for the 
angle that at Q = it holds 

cos(tf) = ^ = /t (4.3) 
\b\ \r\ 



so that 



X max >0=^r>0=^ COS(^) = 1 

$ = 27rm, meZ, 
X m in < r < =► cos(i?) = -1 

=>■ i? = 7T + 27rn, neZ. (4.4) 

Now, the fitting of the various solution branches can 
be performed as follows (see Fig. ©). We start with 
the branch Ti := (R + < R < +oo,X mm < X < 
+00). Moving in on this branch from X w +00 we 
are working in the one-solution sector Q > with 



£(Ti;Q) 



,1/2 



1/3 



,1/2 



1/3 



(4.5) 




x min 









Q-0 












P.I 







-V-a/Pr) 



Q>0 Q<0 Q>0 

FIG. 2: The schematic drawing of the real solution 
branches and the matching points Pi, 2, 3, 4 This figure 
shows that points p2,3 (correspondingly, 6 — tt, 27r) and 
points Pi, 4 (correspondingly, = 0, 37r) have different na- 
ture. So, p2,3 and Pi. 4 we shall call branching points and 
monotonic points (in the sense that function R is mono- 
tonic in the vicinity of these points), respectively. 



until we hit Q = at X = X max . At this point 
P\ := (T%,X = X max ) e Ti we have to perform the 
first fitting. Due to r > we may choose 

R(T 1 ;Q = 0) = 2? 4 / 3 = 2|&| 1 / 3 (4.6) 
so that as simplest parameter choice in (|3.16|) we set 



Pi — (fi, X — Xn 



= 0) i-> i? = 0, fc = 0. 



(4.7) 

Hence, the parametrization for (Ti,Q < 0) will be 
given as 



R(T U Q < 0) = 2|6| 1/3 cos(6>/3). 



(4.8) 



For later convenience, we have replaced here the $ 
from the equations (|3. 15|) . Q3.16P by 9. The reason will 
become clear from the subsequent discussion. We note 
that on this Ti— segment we may set '0 = 8. Let us 
further move on Ti until its end at X m i n , where again 
Q = 0. Because there was no other point with Q = 
on this path, the smoothly changing 9 can at this local 
minimum Pi = V\ D T2 = (X = X m i n , R = only 
take one of the values 9 = ±ir. For definiteness we 
choose it as 9(Pi) = 7r. Hence, it holds 



2|&| 1 / 3 cos( 7 r/3) 



\b\ 1/3 

y/-a/(&y) = R+ 



(4.9) 



as it should hold. For convenience, we may param- 
eterize our movement on the cubic curve by sim- 
ply further increasing 9. This gives for moving on 



T2 — {R+ — R — R— 1 X. n 



< X < X n 



from 



the local minimum at Pi to the local maximum at 
^3 = T 2 n T 3 = (X = X max ,R = R_) a further in- 
crease of 9 by 7r up to 9(Ps) = 2tt. Accordingly, we 
find the complete consistency 



2\b\ 1/3 cos(27r/3) 



-\b\ 1/3 

y/-a/(6j) = R- . (4.10) 



G 



By further increasing 9 up to 9 — 37r we reach the 
point Pi = (X = X mm ,Q = 0)er 3 with 



R(Pi) = 2|6| 1 / 3 cos(3 7 r/3) = -2\b\^ 3 = -2| 



1 1/3 



(4.11) 

Because of r < we can smoothly fit it to the one- 
solution branch 



R(T 4 ,Q)= r + Q 



,1/2 



1/3 



r-Q 



1/2 



1/3 



by setting trivially 

R(P A ) = 2(-\r\) 1 / 3 = 



-2|r 



1/3 



(4.12) 



(4.13) 



Summarizing, we arrived at a very simple and trans- 
parent branch fitting picture, where all the move- 
ment in the three-solution sector can be parameter- 
ized by choosing the effective angle as 9 £ [0, 3tt]. Fi- 
nally, we have to fit this picture in terms of smoothly 
varying 9 € [0, 37r] with the three-solutions Rk from 
(I3.16P . For this purpose we note that the single value 
*& G [0, 7r] in (|3.16p is a projection of our smoothly 
varying 9 E [0, 37r]. Fixing an arbitrary d one easily 
finds the following correspondences 



9(r u &) = 0, 9(T 2 ,d) 
and hence 



R[9(T 2 ,ti)] 



2\b\ 



= 2|6| 
= 2\b\ 



1/3, 



1/3 



1/3 



2tt 



0, 0(T 3 ,0) = 2tt + 

(4.14) 



R 



(fe=0) 



ft 3 



cos 



cos 







JJ[0(T 3 ,0)] = 2|fe| 1 / 3 cos 



3 

- 2tt 
3 

2tt + 15» 



(fe=-l) 



i? 



(fe=i) 



Hi 



= n 2 



(4.15) 



Analogically, we can obtain rules for fitting procedure 
in the case 7 < 0, a > 0. So, all the fitting mechanism 
is clear now and can be used in further considerations. 



V. DYNAMICS OF THE UNIVERSE AND 
SCALARON 

To study the dynamics of the Universe in our model, 
we assume that the four-dimensional metric g in (12.71) 
is spatially flat Friedmann-Robertson- Walker one: 



g = —dt ® dt + a 2 (t)dx ® dx . 



(5.1) 



Thus, scalar curvatures R and R and scalaron 4> are 
functions of time. Therefore, Eq. (|2.9p for homoge- 
neous field <b reads 



• dU 
3H<j> + -— = , 



(5.2) 



where the Hubble parameter H — a/a and the dotes 
denote the differentiation with respect to time t. Po- 
tential U is defined by Eq. (|2.6I) . Because U depends 



2.5 




FIG. 3: The form of the potential 1)2. 6|) as a multi-valued 
function of X = e A< ^ — 1 in the case D = 4, A4 = 0.1, 
7=1 and a = —1. Points Pi, 2, 3, 4 are denned in Fig. |JSJ. 




X = X„„ cosh(0) 
(9<0 



X = X_ cos(0) 

< e < 3x 



X = X„ ln cosh(6<-3;r) 
0>3ir 



FIG. 4: The schematic drawing of Eq. (|5.3[) in the case 
X min > -1. Here, X mM = ^(1/7) (2|a|/3) 3 , X mm = 
(I/7) (2|a|/3) 3 and 6 max is defined by Eq. ijO]! . 



on i? which is a multi- valued function of <^> (or, equiv- 
alently, of X), the potential U is also a multi- valued 
function of X (see Fig. (J3])) 5 . However, our previous 
analysis _shows that we can avoid this problem making 
X and R single- valued functions of a new field 9 (we 
would remind that we consider the particular case of 
«7 < when a < 0, 7 > 0): 



X{9) 



^(1/7) (2|a|/3) 3 cosh(6>) , 6» < 

y^l/7) (2|a|/3) 3 cos(0) , < 9 < 3tt 

-\J(1/j) (2|a|/3) 3 cosh(<9 - 3tt) , 9 > 3tt 

(5.3) 



5 In spite of the divergency of dR/dX in the branching points 
P2,3, the derivatives dU/dX are finite in these points in ac- 
cordance with Eq. (12 . 14D . Moreover, R and X have the same 
values in branching points for different branches. Therefore, 
the branches arrive at the branching points with the same 
values of dU/dX and Fig. J3J clearly shows it. 



7 



and 



R{9) 



2 v / |g|/(6 7 ) cosh(0/3), 
2 yW(6 7 ) cos(0/3) , 0<i 
-2V|a|/(6 7 ) cosh[(0/3) - tt] , 



6> < ; 
< 3vr ; 
6» > 3tt . 
(5.4) 

The function X = X{0) is schematically given in Fig. 

It is necessary to keep in mind that we consider 
the case /' > — > X > — 1. If we demand that 
> — 1 (in opposite case our graphic X(9) will be 
cut into two disconnected parts) then the parameters 
a and 7 should satisfy the inequality: 



x mm = -vCVt) (2M/3) 3 > 



-1 |a| < -7 



1/3 _ 

The maximal value of 9 (which is greater than 3-7T in 
the case X m i n > —1) is defined from the transcenden- 
tal equation 



(c + Vc 2 - 1 ) V3 + (c - Vc 2 - 1 ) 



1/3' 



cosh[(0 ma2; /3) - tt] 



(5.6) 



where c 



(1/7) (2M/3)' 



— 1 corresponds to the limit 
of Eq. (|5.3p and formula 



The limit X -s- 
x . With the help 



Therefore, the extrema correspond cither to the so- 
lutions of the equation (|2.10[) dU /d(f> = for finite 
dX/dO (X > —1) or to the solutions of the equation 
dX/de = {X > -1) for finite dU/d(j>. The form 
of the potential U ( see Eqs. (|2.4I) and (|2.6l0 shows 
that this potential and its derivative dU/dd is finite 
for X > -1. Thus, as it follows from Eq. (|5.8[) . the 
potential J7(0) has extrema at the matching points 
9 = 0, 7r, 27r, 37r. Additional extremum points are real 
solutions of the equation (|2 . 10|) . For our model (|3 . 1 1) 
this equation reads: 



i? 4 7 



R 2 a 



D 



- 2 



R 



D 



- 1 - DAn = 



(5.11) 



The form of this equation shows that there are two 
particular cases: D = 8 and D = 4. The D = 8 case 
was considered in |ll|. Let us consider now the case 
D = 4: 



i? 4 



2 7 7 



0. 



(5.12) 



It is worth of noting that parameter a disappeared 
from this equation. Thus a has no an influence on 
a number of additional extremum points. To solve 
this quartic equation, we should consider an auxiliary 
cubic equation 



dX 



d9 A(X + 1) d9 
we can also get the following useful expressions: 



(5.7) 



d9 



dX 

~dl 



= 0. 



(5.8) 



= 0,7T,27r,37T 



It can be easily verified that field 
equation: 



satisfies the 



9 + 3H9 + T{9)9 



52 



(5.9) 



Here, we introduce the one dimensional metric on the 
moduli space G(9) = G 11 = (Gn)" 1 = {d4>/d9)- 2 
with the corresponding Christoffel symbol T(9) = 
T\ t = (l/2)G 11 (G 11 ), e = (d 2 ^/d9 2 )/(d^/d9). 



8A4 1 

u - = ■ 

7 47^ 



(5.13) 



The analysis of this equation can be performed in sim- 
ilar manner as we did it for the cubic equation (I3.3[) . 
Let us introduce the notations: 



A4 
3 7 



1 1 



Q:-- 



-3 1 

9 =7 



8A4 



.(5.14) 



It make sense to consider two separate cases. 

1. sign 7 = — signA 4 => Q > . 

In this case we have only one real solution of Eq. 
(EH: 



Ml 



jl/2 



1/3 



/ 1 1/3 

Q 1/2 >0. 



(5.15) 



A. Properties of the potential U(8) 



Then, solutions of the quartic (I5.12p are the real roots 
of two quadratic equations 



As we mentioned above, the potential (|2.6|) as a 
function of 9 is a single- valued one. Now, we want 
to investigate analytically some general properties of 
U{9). In this subsection, D is an arbitrary number of 
dimensions and signs of a and 7 are not fixed if it is 
not specified particularly. 

First, we consider the extrema of the potential U{9). 
To find the extremum points we solve the equation 



dU 

~dJ 



dU d<t> 
l6d9 



dU 



dX 



A(X + 1) d9 



0. 



(5.10) 



R 2 ± y/uiR+ i Fiti ± e-\J u\ + 3qJ — , e = sign7. 

(5.16) 

Simple analysis shows that for any sign of 7 we obtain 
two real solutions: 



7 < =*> R\ 



7 > r[J 



(+) _ 1,1/2 



l /2 ± v /-in 1 + i(nf + 3g)V2 : 



1 1/2 

2 U 



^i^-i^ + iK + Sg)^. 

(5.17) 
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2. sign 7 = signA4 =>■ Q ^ . 
It is not difficult to show that in this case the real 
solutions of the form of (|5.17l) (where we should make 



u 



the evident substitution 
place if 



Q >0 



hi 



1/3 



3q — > u\ — i\q\) takes 
<^J^- (5.18) 



32|A 4 



Now, we want to investigate zeros of the potential 
U. For /' ^ => X ^ -1, the condition of zeros of 
the potential (I2.4[) is: 

#/'-/ = => + aR 2 + 2A D = 0. (5.19) 

Therefore, zeros are defined by equation: 



| 

67 



2Ap 

37 



-1 1/2 



(5.20) 



Obviously, the necessary conditions for zeros are: 
7>0 =*> A D < a 2 /(24 7 ) , 

7<0 A D > -a 2 /(24| 7 |). 



(5.21) 



Additionally, we should check that the r.h.s. of the 
equation (|5.20l) is positive. 

Let us consider now asymptotical behavior of the 
potential U(9). Here, we want to investigate limits 
9 — > 9 m ax and 9 — > — 00. In the former case we get: 



9^9 max U{6) -> -sign (f(6 r , 

In the latter case we obtain: 



„)) x 00. 
(5.22) 



->• -00 =>• 



U{9) ~ exp 



8- £> 
D-2 



+00 , D > 8 
const > , D = 8 
+0 , 2 < £> < 8 



(5.23) 

where we used Eqs. (|5.3[) and (|5.4j) . This result coin- 
cides with conclusions of Appendix A in [16[ . 

To illustrate the described above properties, we 
draw the potential U(9) in Fig. §5§ for the follow- 
ing parameters: D — 4, A4 = 0.1, 7 = 1 and a = — 1. 
These parameters contradict to the inequalities (|5.18p 
and (|5.2ip . Therefore, 9 = 0,7r, 2tt, 3tt are the only 
cxtremum points of the potential U (9) and zeros are 
abse nt. T hese parameters are also satisfy the condi- 
tion (|5.5I) . The absence of zeros means that all minima 
of the potential U{9) are positive. 

For our subsequent investigations, it is useful also 
to consider an effective force and mass squared of the 
field 9. As it follows from Eq. (|5.9p . the effective force: 



(5.24) 



Varying Eq. (|5.9[) with respect to field 9 we obtain dy- 
namical equation for small fluctuation 89 where mass 
squared reads: 



G(9) 



d 2 U dG{9) dU 
~d!P + d9 ~d9 



(5.25) 



2 . 5 




-2 71 



2 7T 



3 n 



G 



FIG. 5: The form of the potential (|2.6|l as a function of 
9 in the case D — 4, A4 = 0.1, 7=1 and a = — 1. For 
these values of the parameters, all extrema correspond to 
the matching points 8 — 0, it, 2n, 37r. In the branching 
points 9 = 7r, 2tt the potential has local maximum and 
local non-zero minimum, respectively, and the monotonic 
points 9 — 0, 3tt are the inflection ones. Potential tends 
asymptotically to +00 when 9 goes to 9 m ax and to zero 
when 9 — s> — 00. 



In Fig. (J6j> we show the effective force and the mass 
squared as functions of 9 for the potential drawn in 
Fig. ([5]). These figures indicate that field 9 may 
have very nontrivial behavior. This non-triviality fol- 
lows from two reasons. First, the field 9 has non- 
canonical kinetic term which result in appearing of 
non-flat moduli space metric G(9) and derivative of 
G(9) in Eq. Second, the function G{9) has sin- 

gular behavior at the matching points 9 = 0, 7T, 27r, 37T. 
Thus, our intuition does not work when we want to 
predict dynamical behavior of fields with equations of 
the form of (|5.9[) with potential drawn in Fig. ([5]), 
especially when field approaches the matching points. 
It is necessary to solve equations analytically or to 
investigate them numerically. Such analysis for our 
model is performed in the next subsection where we 
concentrate our attention to the most interesting case 
where all extrema correspond to the matching points 

9 = 0,7T,27T,37r. 



B. 



Dynamical behavior of the Universe 
and field 9 



Now, we intend to investigate dynamical behavior 
of scalar field 9 and the scale factor a in more de- 
tail. There are no analytic solutions for considered 
model. So, we use numerical calculations. To do it, 
we apply a Mathematica package proposed in [26| and 
adjusted to our models and notations in Appendix A 
of our paper [l8j]. According to these notations, all di- 
mensional quantities in our graphics are given in the 
Planck units. Additionally, in the present paper we 
should remember that metric on the moduli space is 
not flat and defined in Eq. (|5.9p . For example, the 
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FIG. 6: The effective force (|5.24[1 (left panel) and the mass squared (|5.25[) (right panel) for the potential U(ff) drawn in 
Fig. ([5]). These pictures clearly show singular behavior of F and nig in the matching points 6 = 0, it, 2tt, 3n. 



canonical momenta and the kinetic energy read: 

where 8ttG = k\ and G is four-dimensional Newton 
constant. To understand the dynamics of the Uni- 
verse, we shall also draw the Hubble parameter: 

3^ = 'SH 2 = \GiiB 1 + U{6) (5.27) 
and the acceleration parameter: 

Fig. © shows that the effective force changes its 
sign and the mass squared preserves the sign when 9 
crosses the branching points it, 2tt and vise verse, the 
effective force preserves the sign and the mass squared 
changes the sign when 9 crosses the monotonic points 
0, 3-7T. Therefore, it make sense to consider these cases 
separately. 



1. branching points 8 — tt,2tt 

First, we consider the dynamical behavior of the 
Universe and scalaron in the vicinity of the branch- 
ing point 9 = 2ir which is the local minimum of the 
potential in Fig. ([5]). The time evolution of scalaron 
field 9 and its kinetic energy Eki n are drawn in Fig. 
([7]). Here and in all pictures below, we use the same 
parameters as in Fig. ([3]). The time t is measured 
in the Planck times and classical evolution starts at 
t = 1. For the initial value of 9 we take 9i n iu a i = 3.5. 
We plot in Fig. [5] the evolution of the logarithms 
of the scale factor a(t) (left panel) and the evolution 
of the Hubble parameter H(t) (right panel) and in 
Fig. IH1 the evolution of the parameter of acceleration 



q(t) (left panel) and the equation of state parameter 
u)(t) = [2q(t) + l]/3 (right panel). 

Fig. ([7]) demonstrates that scalar field 9 bounces an 
infinite number of times with decreasing amplitude in 
the vicinity of the branching point 9 — 2tt. 9 cannot 
cross this point. From Figs. © and © we see that 
the Universe asymptotically approaches the de Sitter 
stage: H — s- const, q — s- +1, and u — > —1. Such 
accelerating behavior we call bouncing inflation. 

Concerning the dynamical behavior in the vicinity 
of the branching point 9 = it, our analysis (similar per- 
formed above) shows that the scalaron field 9 cannot 
cross this local maximum regardless of the magnitude 
of initial velocity in the direction of 9 — it. It bounces 
back from this point. 



2. monotonic points 9 = 0, 3tt 

Now, we want to investigate the dynamical behavior 
of the model in the vicinity of the monotonic points 
9 = 0, 37r which are the points of inflection of the po- 
tential in Fig. ©. Figs. ([5]) and ^ show that for 
both of these points the model has the similar dy- 
namical behavior. Therefore, for definiteness, we con- 
sider the point 9 = 3tt. To investigate numerically the 
crossing of the monotonic point 3-7r, it is necessary to 
take very small value of a step At. It can be achieved if 
we choose very large value of the maximum number of 
steps. Thus, for the given value of the maximum num- 
ber of steps, the closer to 37r the initial value 9i n iu a i 
is taken the smaller step At we obtain. For our cal- 
culation we choose 9 in iti a i = 9.513. Fig. (|TU)) demon- 
strates that scalar field 9 slowly crosses the monotonic 
point 37T with nearly zero kinetic energy 6 . Then, just 
after the crossing, the kinetic energy has its maximum 
value and starts to decrees gradually when 9 moves to 
the direction 2tt. Figs, (jlll) and (fT2j) demonstrate 
the behavior of the Universe before and after crossing 



The derivative dd/dt goes to —00 when 8 — > 3tt (with different 
speed on different sides of 3n) but d<p/d9 = at 37r and kinetic 
energy is finite (see Eq. J5J26}). 
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FIG. 7: Dynamical behavior of scalar field 9(t) (left panel) and its kinetic energy Ekmit) (right panel) in the vicinity of 
the branching point 8 = 2n. 
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FIG. 8: The time evolution of the logarithms of the scale factor a(t) (left panel) and the Hubble parameter H (t) (right 
panel) for the trajectory drawn in Fig. ([7]). A slightly visible oscillations of ln(a) (caused by bounces) can be seen by 
magnification of this picture (middle panel). 



37r. We do not show here the vicinity of the branching 
point 2ir because when 9 approaches 2tt the Universe 
has the bouncing inflation described above. Hence, 
there are 3 phases sequentially: the short de Sitter- 
like stage during slow rolling in the vicinity of the 
inflection point before crossing, then decelerating ex- 
pansion just after the crossing with gradual transition 
to the accelerating stage again when 9 approaches the 
branching point 2n. Clearly that for another mono- 
tonic point 9 = we get the similar crossing behavior 
(without the bouncing stage when 9 — > — oo). There- 
fore, the monotonic points 9 = and 9 = 37r are 
penetrable for the scalaron field 9. 



VI. CONCLUSIONS 

We have investigated the dynamical behavior of the 
scalaron field <j> and the Universe in nonlinear model 
with curvature-scmared_and curvature-quartic correc- 
tion terms: f(R) = R + aR 2 + 7 i? 4 - 2A D . We 
have chosen parameters a and 7 in such a way that 
the scalaron potential U(<fr) is a multi- valued function 
consisting of a number of branches. These branches 
arc fitted with each other either in the branching 
points (points i-2,3 in Fig- ©) or in the monotonic 
points (points in Fig. ([3])). We have reparam- 
eterized the potential U in such a way that it be- 



comes the one-valued function of a new field vari- 
able 9 = 9((f>) (see Fig. ([5])). This has enabled us 
to consider the dynamical behavior of the system in 
the vicinity of the branching and monotonic points in 
(D = 4)-dimensional space-time. Our investigations 
show that the monotonic points are penetrable for 
scalaron field (see Figs. (TTUl) - (fl"2)) ) while in the vicinity 
of the branching points scalaron has the bouncing be- 
havior and cannot cross these points. Moreover, there 
are branching points where scalaron bounces an infi- 
nite number of times with decreasing amplitude and 
the Universe asymptotically approaches the de Sitter 
stage (see Figs. ©-©). Such accelerating behavior 
we call bouncing inflation. It should be noted that for 
this type of inflation there is no need for original po- 
tential U {4>) to have a minimum or to check the slow- 
roll conditions. A necessary condition is the existence 
of the branching points. This is a new type of infla- 
tion. We show that this inflation takes place both in 
the Einstein and Brans-Dicke frames. We have found 
this type of inflation for the model with the curvature- 
squared and curvature-quartic correction terms which 
play an important role during the early stages of the 
Universe evolution. However, the branching points 
take also place in models with i? _1 -type correction 
terms [19) . These terms play an important role at 
late times of the evolution of the Universe. There- 
fore, bouncing inflation may be responsible for the 
late-time accelerating expansion of the Universe. 
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FIG. 9: The parameter of acceleration q(t) (left panel) and the equation of state parameter u(t) (right panel) for the 
scale factor in Fig. ((8|. 



e de/dt Ek 




FIG. 10: Dynamical behavior of scalar field 8(t) (left panel) and its time derivative d9/dt (middle panel) and kinetic 
energy Eki n (t) (right panel) for the case of crossing of the inflection point — 3tv. 



To conclude our paper, we want to make a few com- 
ments. First, there is no need for fine tuning of the ini- 
tial conditions to get the bouncing inflation. In Figs. 

we have chosen for definiteness the initial con- 
ditions 9 = 3.5 and Eki n = 0. However, our calcu- 
lations show that these figures do not qualitatively 
change if we take arbitrary 9 G (ir, 2ir) and non-zero 
Ekin- Second, Figs. ^ indicates that the minimum 
at 9 = 2ir is stable with respect to tunneling through 
the barrier at this point. The situation is similar to 
the quantum mechanical problem with infinitely high 
barrier. We have already stressed that the form of 
the potential Fig. §5§ is not sufficient to predict the 
dynamical behavior of 9. This field has very non- 
trivial behavior because of non-canonical kinetic term 
and singular (at the matching points) non-flat mod- 
uli space metric G(9). Therefore, it is impossible to 
"jump" quantum mechanically from one branch to an- 
other. We cannot apply to our dynamical system the 
standard tunneling approach (e.g., in [27|). This prob- 
lem needs a separate investigation. Third, it is worth 
noting that the Universe with a boun ce p receding the 
inflationary period was considered in |28| where it was 
shown that due to a bounce the spectrum of primor- 
dial perturbations has the characteristic features. It 
indicates that the similar effect can take place in our 
model. This is an interesting problem for future re- 
search. 
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Appendix A: Bouncing inflation in the 
Brans-Dicke frame 

According to Eq. (|2.7I) . the four-dimensional FRW 
metrics in the Einstein frame (|5.ip and in the Brans- 
Dicke frame are related as follows: 

-dt®dt+a 2 \t)dx®dx = f [—dt <g> dt + d 2 (t)dx <g> dx] , 

(A.1) 

wher e / = X + 1 > and X is parameterized by Eq. 
(|5.3I) . Therefore, for the synchronous times and scale 
factors in both frames we obtain, correspondingly: 

dt = dt/y/f 7 ®, (A.2) 
d(t) = a{t{t))l^nW), (A.3) 
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FIG. 11: The time evolution of the loga rith ms of the scale factor a(t) (left panel) and the Hubble parameter H(t) (right 
panel) for the trajectory drawn in Fig. (|10|l . 
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FIG. 12: The para mete r of acceleration q(t) (left panel) and the equation of state parameter u>(t) (right panel) for the 
scale factor in Fig. (1111) . 
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FIG. 13: The time evolution of the logarithms of the scale factor a(t) (left panel) and the Hubble parameter H(t) (right 
panel) for the trajectory drawn in Fig. (J7J) in the Brans-Dicke frame. 



which lead to the following equations: 



dt 



where we choose the constant of integration in such a 
way that i(t = 1) = 1, and 



H{t) 



da 1 
dt d 



V* (*(*)) + ! 



H(t(t}) 



dX 



2(X(t{t)) + 1) dt 



(*(*)) 



(A.5) 



(A 4) From the latter equation we get the relation between 
the Hubble parameters in both frames. We plot in 
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Fig. (fT3"]) the logarithms of the scale factor a(t) and 
the Hubble parameter H(t) for the trajectory drawn 
in Fig. (J7J . These pictures clearly demonstrate that in 
the Brans-Dicke frame the Universe has also asymp- 
totical de Sitter stage when the scalaron held ap- 
proaches the branching point 9 = 2n. It is not dif- 
ficult to verify that because X(t — > +00) — > X max 



and dX/dt(t — > +00) — > 0, we obtain the following 
relation for the asymptotic values of the Hubble pa- 
rameters in both frames: 



H = Hy/X max + 1 . (A.6) 
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